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Abstract: We introduce a new method for applying an inflationary prior to a cosmological 
dataset that includes relations between observables at arbitrary order in the slow roll 
expansion. The process is based on the inflationary flow equations, and the slow roll 
parameters appear explicitly in the cosmological parameter set. We contrast our method 
to other ways of imposing an inflationary prior on a cosmological dataset, and argue that 
this method is ideal for use with heterogeneous datasets. In particular, it would be well 
suited to exploiting any direct detection of fundamental tensor modes by a BBO-style 
mission. To demonstrate the practical use of this method we apply it to the WMAPI+All 
dataset, and the newly released WMAPII dataset on its own and together with the SDSS 
data. We find that all basic classes of single field inflationary models are still allowed at 
the 1 — 2(7 level, but the overall parameter space is sharply constrained. In particular, we 
find evidence that the combination of WMAPII+SDSS is sensitive to effects arising from 
terms that are quadratic in the two leading-order slow roll parameters. 
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1. Introduction 

Observational cosmology is a booming science. We typically see at least one significant 
data release each year, from both Cosmic Microwave Background [CMB] experiments and 
probes of Large Scale Structure [LSS]. Several generations of experiments are now being 
designed and constructed, suggesting that current rates of progress will continue well into 
the future. A major and very welcome challenge for the cosmological community is how 
to best exploit this torrent of data. Each major data release is typically accompanied by 
newly estimated values of the standard cosmological parameter set - based either on the 
new dataset alone, or an updated "world average" drawn from the combination of several, 
complementary experiments. 

There are two ways to construct a cosmological parameter set - the first and more 
general is to use a purely empirical set of observables which contains as few theoretical 
assumptions as possible. Alternatively, one may impose a "prior" corresponding to the 
predictions of a specific theoretical model. Because this reduces the overall parameter 
space, this class of variables can be more tightly constrained than a generic set. The 
most common theoretical prior is derived from inflation. This typically puts two cuts 
on the parameter set. Firstly, without careful tuning, inflation predicts that present day 
universe is flat, so we impose ilTotai = 1- Secondly, generic models of single field, slow 
roll inflation obey a consistency condition that relates the amplitudes of the primordial 
scalar and tensor spectra with the slope of the tensor spectrum (see Q for a discussion.) 
In addition to these two constraints, the conventional description of the primordial spectra 



- 1 - 



in terms of their amplitude A and spectral index n is unnecessary, if one assumes that the 
primordial perturbations are inflationary.^ 

The spectral indices are derived from a truncated Taylor expansion of the overall power- 
spectrum, P{k), where k is the comoving wavenumber. One may extend this expansion to 
higher order by including the running, a = dng/dlnk, as an additional parameter. This 
is obviously an empirical characterization of the spectrum. Inflationary theory expresses 
these parameters as functions of the potential V{4>) and its derivatives. This introduces 
a further truncated expansion ~ this time in terms of the slow roll parameters, which are 
dimensionless combinations of derivatives of the potential. This expansion is independent of 
the Taylor expansion underlying the spectral indices, and ignores the relationships between 
the parameters of the spectra found in the slow roll limit. Consequently, once we have 
adopted an inflationary prior there is nothing to be gained by expressing the spectrum 
in terms of Us and dug/dlnk, and doing so necessarily invokes a further and ultimately 
unnecessary truncated expansion and risks discarding useful information. In this paper we 
will investigate the use of a set of cosmological parameters which incorporates the slow roll 
variables directly, making no direct use of the spectral index or its running. 

If all we gained from this substitution was a reshuffling of the variables that describe 
the perturbation spectrum, the issue would be of purely academic interest. However, 
the slow roll parameters have an immediate physical meaning in inflationary models, and 
including them in our parameter set raises a number of new and exciting possibilities. In 
conventional fits to the spectrum, one specifies the amplitude and index at some fiducial 
wavenumber, ko (0.002Mpc~^ in WMAP data), and the slow roll parameters are likewise 
measured at some fixed wavenumber. The simplest way to implement this in a CMB code 
such as CAMB or CMBFast is to express and dus/dlnk as a function of the slow roll 
parameters at ko, and then compute the CMB and LSS power spectra in the usual way. 
However, the slow roll parameters are not constant throughout the inflationary epoch but 
are fully specified by the inflationary flow equations a set of coupled flrst order 

differential equations. The values of the slow roll parameters at /cq provide the initial data 
for this system, and solving the flow equations yields values of the slow roll parameters 
at any other value of k. Formally, the slow roll hierarchy is inflnite, but any dataset will 
only contain enough information to put meaningful constraints on a finite subset of terms. 
However, if we truncate the slow roll hierarchy at any given point, the structure of the 
equations ensures that the truncation holds everywhere, and does not depend explicitly on 
the fiducial wavenumber kg. 

The flow equations contain enough information for us to reconstruct the inflationary 
potential up to an overall normalization - which would only be fixed by an unambiguous 
measurement of the tensor contribution to the CMB. Consequently, this approach provides 
the most direct route to a long-held ambition of early universe cosmology: recovering the 
inflationary potential directly from astrophysical data. Previous attempts to do this have 

^As opposed to both "defect" models such as cosmic strings, and also non-standard inflationary schemes, 
such as the curvaton scenario [^. In the latter case, inflation solves the "global" cosmological problems 
such as homogeneity, isotropy and flatness, but the perturbation spectrum is not immediately determined 
by the inflationary potential. 
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typically relied on computing the spectral parameters, inverting these to obtain the slow 
roll parameters that gave rise to them, and then solving these to obtain the underlying 
potential. This process can be carried out analytically, as summarized in [p. To first order, 
the mappings between spectral indices and the potential slow roll parameters are relatively 
simple, but they become considerably more complicated at second order By working 

directly with the Hubble slow roll hierarchy, our algorithm automatically incorporates 
any and all relations between the inflationary observables, to whatever order we choose 
to truncate the expansion. In particular, we will see that the new WMAP dataset, in 
combination with SDSS data, has enough power to significantly constrain the allowed 
running of the scalar spectral index. This arises from terms proportional to the square 
of the first two slow roll parameters, and is thus indicative of both the quality of present 
cosmological data, and the ability of our method to exploit this by including higher order 
relations between the slow roll parameters and inflationary observables. 

There are two important precedents for this work. The first is Monte Carlo Recon- 
struction and the second is Leach et aZ.'s progress towards placing direct constraints of 
the slow roll parameters using both CMB and LSS data |^, |^, |lO[. As initially proposed, 
Monte Carlo reconstruction proceeds by solving the flow equations with arbitrary initial 
conditions and then identifying models which match some predefined constraint (s) on the 
spectral parameters. The problems with this approach is that existing implementations a) 
consider a vast number of models which do not provide a good fit to data, b) there is no 
natural measure with which to weight the choice of initial values of the slow roll parameters 
and c) no attempt to establish a relative likelihood for the parameter values that do satisfy 
the chosen astrophysical constraints. All these problems are solved when we obtain the 
slow roll parameters directly from a Monte Carlo Markov Chain [MCMC] flt to data. One 
of the virtues of the MCMC approach is that one can draw parameter values from a flat 
distribution, and the adaptive nature of the chains ensures that you focus computational 
effort on the subspace where the likelihood is most sharply peaked. Moreover, the relative 
merit of different parameter choices is quantifled by the likelihood computed at each step 
in the chain with respect to the data. 



Conversely, Leach et al. fg, y, |10[, work with rig and dug/dlnk, as first suggested 



in [^]. The cosmological observables are expressed in terms of the slow roll parameters, 
whose scale dependence is included via the slow roll approximation. While this approach 
could, in principle, be extended to arbitrary order in slow roll, the algebraic complexity 
encountered at each successive order is significantly higher than the last ||l^ . On the other 
hand, because only a comparatively small piece of the inflationary potential is directly 
"sampled" by observations of the CMB and LSS, in most cases the second order results 
are sufficiently accurate for practical purposes. Indeed, in what follows, we compare our 
approach to that of Leach et al. when applied to the WMAPI+All dataset, and show that 
it leads to essentially identical constraints on the slow roll parameters. On the other hand, 
the inevitable improvement of cosmological data will increase both the range of wavelengths 
over which we can probe the spectrum, and the accuracy to which we can measure it. If one 
ever wishes to exceed the limitations of the second order expansion, then the flow equations 
provide an elegant route to reconstruction. In particular, if primordial gravitational waves 
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are ever detected directly via a BBO-style experiment the range of scales over which we 
can access primordial perturbations would increase dramatically, and the flow equations 
would be immediately applicable to this scenario (see ||l^, 14, 15, |l6|).^ In the absence 



of a direct detection of primordial gravitational waves at solar system scales, future data 
will still extend the range of k values over which we can probe the primordial spectrum - 
in particular, close quasar pairs facilitate short scale measurements of the Lyman-a forest 
power spectrum [O], and future measurements of the high redshift 21 centimeter radiation 



will probe the primordial power spectrum at scales up to /c ~ lOMpc"^ [18|. 

A final advantage of this approach is that it allows us to impose constraints based on 
the number of e-folds that elapse between some fiducial mode with wavenumber crossing 
the horizon and the end of inflation. In their work, Leach et al. note that because they 
extrapolate the power spectrum about a point there is a possibility that it will cross zero in 
the range of interest, leading to unphysical values for the power spectrum. In our approach, 
we see that for some parameter values e becomes larger than unity in the region of interest, 
signaling the end of inflation. In this case we can add a prior that inflation lasts longer 
than some minimal number of e-folds and exclude the model on physical grounds. 

One of the virtues of the Hubble slow roll formalism is that, if it is used self consistently, 
the resulting inflationary dynamics actually provide an analytic solution to the evolution 
equations for a scalar field dominated FRW universe. However, the perturbation spectrum 
can typically only be obtained exactly in very special circumstances |19, 20, 21], so even with 
an exact result for the background one typically has to use the slow roll approximation in 
order to obtain the perturbation spectra. Grivell and Liddle describe a reconstruction 
technique that again avoids the use of the spectral indices when computing inflationary 
perturbation spectra, and goes one step further in that it entirely abandons the use of the 
slow roll approximation, directly integrating the fundamental mode equations in order to 
obtain the power spectrum. Here we employ slow roll expressions for the perturbation 
spectrum itself, accurate to second order in slow roll. This may appear paradoxical, but 
since we are accounting for the scale dependence in the slow roll parameters by using the 
flow equations, the higher order corrections to the potential are fed into Pr and Ph via the 
scale dependence of e and ry. However, it would be entirely feasible to evaluate the spectrum 
exactly (up to the assumptions inherent in the use linear gravitational perturbation theory) 
in our approach. 

Finally, we believe that using the HSR parameters directly in the chains avoids the 
theoretical uncertainties in the inflationary predictions associated with a given potential, 
discussed by Kinney and Riotto |^]. These arise from the ambiguity in the matching 
between present physical scales and the number of e-folds before the end of inflation at 
which that scale left the horizon. Since we are expanding the potential around some fiducial 
point, we never need to explicitly calculate the matching, and solve an exact relationship 
for the connection between In A; and the inflaton field (p. Furthermore, we can check that 
any given potential permits some minimum number of e-folds to elapse after cosmological 
scales have left the horizon. On the other hand, any such calculation is likely to be inexact 

^BBO would detect inflationary perturbations at scales of a few million kilometers, rather than mega- 
parsecs. 
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since the higher slow roll terms can easily be significant at the end of inflation - and 
inflation may end with a hybrid style transition induced by an instability in a direction in 
field space that is orthogonal to the inflation. 

The key assumption that underlies this work is that the primordial perturbations were 
generated by the standard inflationary mechanism - and that during the epoch in which 
the perturbations are generated, the evolution of the universe can be accurately modeled in 
terms of a single, slowly rolling, minimally coupled scalar field. We expect that if there is 
an accurate description of the early universe of this form, it will be via an effective theory, 
rather than one arising from the explicit presence of an appropriate field in the Lagrangian 
that describes GUT scale physics. In particular, there is a relation between the overall 
excursion made by the inflaton field and the energy scale of inflation [^, |2^. This can be 
used to argue, on the basis of effective field theory considerations, that it is very difficult 
for the inflaton to be a single, fundamental scalar field embedded within string theory. 
However, models relying on multiple scalar fields do arise naturally in string theory (see 
e.g. ||2^, ^), and the method described here will cover the effective single field description 
of these models, assuming that one exists. Lastly, by using a truncated slow-roll hierarchy, 
we are implicitly assuming that the inflationary potential has no sharp "features" in the 
region traversed by the inflaton as observationally relevant scales leave the horizon. Models 
of this type were explicitly compared to CMB data in [^], and the computation of the 
perturbation spectrum in this situation frequently exceeds the capacity of the slow roll 



approximation |29] 



This paper serves two purposes. On the one hand, it explores a new approach to 
extracting the most leverage from applying an inflationary prior to a cosmological dataset, 
while it also addresses the more topical question of the inflationary constraints that can be 
gleaned from the WMAPII dataset. We present chains for the WMAPI+All dataset, which 
is comprised of the original WMAPI release 0], ACBAR ||, CBI ||, VSA ||, the 



HST prior on the Hubble constant |35|, and a tophat age prior that the universe is between 
10 and 20 Gyr old. We use the WMAPI+All chains to compare the different parameter 
estimation methodologies, and to verify our codes. We then apply the HSR reconstruction 
algorithm to the WMAP dataset ||, |8| on its own and to WMAPII+SDSS.^ We will 
find that for these datasets, the second order expressions for the running of the slow roll 
parameters (as in ||8|, |^, yield essentially identical results to those obtained the full 
Hubble slow roll formalism. However, this work shows that the flow equations can be used 
directly in a Monte Carlo Markov Chain based parameter estimation and, as noted above, 
the HSR approach is ideally to the high quality, heterogeneous datasets we expect to be 
dealing with in the future. Our specific results confirm that the combination of WMAPII 
and SDSS puts strong constraints on the infiationary parameter space, whereas WMAPII 
on its own is roughly equivalent to the WMAPI+All dataset. 

The rest of this paper is organized as follows. In the next section, we review the 
Hubble slow roll formalism and the computation of the perturbation spectra. In Section ^ 
we review the convergence issues and assumptions implicit in the use of a truncated slow 



^Sloan Digital Sky Survey M 
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roll hierarchy, in Section Q we describe the implementation of a Markov Chain Monte Carlo 
analysis with this parameter set, and give our results in Section ^ for both the WMAPI+All 
dataset, and the new WMAPII data |36, 37, 38 1. We discuss our constraints on the slow 
roll parameters in Section ^ and conclude in Section |^ 



2. Slow Roll Hierarchy and the Inflationary Potential 

The inflationary flow hierarchy provides a convenient method for generating potentials 
satisfying slow roll conditions for use with Monte Carlo techniques. In this picture, one 
takes the view that the evolution of spacetime is governed by a single order parameter, 
which can be written in form of a minimally coupled scalar field. Under this assumption, 
one can reformulate the exact dynamical equations for inflation as an infinite hierarchy 
of flow equations described by the gener ahzed "Hubble Slow Roll" (HSR) parameters g 
^, 0, |2^, 39]. In the Hamilton- Jacobi formulation of inflationary dynamics, one expresses 
the Hubble parameter directly as a function of the field (p rather than a function of time, 
H = H[(j)), under the assumption that (j) is monotonic in time. The equations of motion 
for the field and background are given by: 
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(2.2) 



PI 



Overdots correspond to time derivatives and primes denote derivatives with respect to (p. 
Equation referred to as the Hamilton- Jacobi Equation, allows us to consider inflation 
in terms of H{(j)) rather than V{4>). The former, being a geometric quantity, describes 
inflation more naturally. 

The major advantage of this picture is that it removes the field from the dynamics, 
and lets one study the generic behavior of slow roll inflation without making detailed 
assumptions about the underlying particle physics. The HSR parameters ^Xh are defined 
by the infinite hierarchy of differential equations 
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Substituting Eq. 2.3 into Eq. immediately gives the potential in terms of e: 
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The flow equations allow us to consider the model space spanned by inflation. The 
trajectories of the flow parameters are governed by a set of coupled first order differential 
equations. As first pointed out by Liddle |4C], these equations have an analytic solution. 
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Truncating the hierarchy of flow parameters, so that the last non-zero term is Xh ensures 
that *^+^A// = at ah times. From Eq. O it fohows that 







(2.6) 



at all times. We thus arrive at the general solution H((j)) as an order M + 1 polynomial in 

M+l 



Hi 



l + Bi 



mpi 



+ 



B 



M+l 



mpi 



(2.7) 



So far, we have not specified an initial value for cj). If we truncate the series at ^ Xh-, 
any set of flow parameters spans an M + 1 dimensional space. However, the evolution 
equations evolve the ^ Xh backwards and forwards as a function of Consequently, the 
set of distinct trajectories spans an M dimensional space, effectively flbering the space of 
initial conditions for the flow hierarchy. However, if the flow parameters are specified at 
(/) = the ambiguity is removed. 

Further, from the definition of e((/)), 



m 



PI 



47r 



\mpi 



+ --- + (M + 1) 



M+l \ / </> 



mpi 



\ mpi 



M 



^ \ mpi 



H + Bm+1 



mpi 



Af+1 



The coefficients Bi can be written in terms of the initial values of the HSR parameters as 
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The sign of Bi specifies in which direction the field is rolling. Without loss of generality, 
we can pick some fiducial physical scale that corresponds to (/> = (which we choose to be 
/cq = 0.002 Mpc~^). Then, with the above convention, > corresponds to scales larger 
than /cq, and < corresponds to smaller scales. We associate a physical wavenumber 
with a value of by solving the additional differential equation, 
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Finally, we must express the scalar and tensor primordial power spectra in terms of these 
parameters. A first order expansion around an exact solution for the case of power law 
inflation ll| gives 
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where r] = ^Xh, C = -2 + ln2 + 7« -0.729637 and 7 is the Euler-Mascheroni constant. 
One minor inconsistency in our approach is that ( 2.12| ) and (2.13) are calculated using the 



slow roll approximation, and truncated at second order. In practice, we do not believe 
this is a significant drawback, as the datasets we are working with provide meaningful 
constraints on (at most) the lowest three slow-roll parameters. In the future, this problem 
could be surmounted by either using a higher order expansion for the spectra, or even 
integrating the mode equations directly. 

Now we have all the ingredients to use this formalism for reconstructing the inflationary 
potential, using a Markov Chain Monte Carlo (MCMC) analysis. Once we choose the order 
M at which to truncate the HSR series, we pick the values of the HSR parameters at the 
fiducial scale, giving H and all the HSR parameters as functions of (j). After solving for 
A;((/>) we can map our results between k and (j), giving us the primordial power spectra as a 
function of k, from which we can then compute the CMB power spectra and estimate the 
cosmological parameters in the usual way. Finally, after constraining the HSR parameter 



values via a MCMC analysis, Eq. |2.5| can be used to reconstruct the corresponding family 
of potentials. Finally, recall that H{(j)) has an overall multiplicative normalization, Hq. 
Within the slow roll hierarchy this is a free parameter. However, once we specify the Bi at 
ko, Hq is fixed by the amplitude of the observed scalar perturbation spectrum, Ag, after 
inverting equation p.l2 .^ 



3. Convergence Issues and Implicit Assumptions 

It is worth noting is that we are effectively fitting to a specific functional form for the 
potential: once we specify -ff(</>) we have also specified V{(j)). To see this more clearly, let 
us cut off the slow roll expansion at B2. In this case, 

2 tt2 

Vi(j)) = ^^4^ r(127r - Bf) + 4(6^^1 - BiB2)f + 4(3^Bi + 671^2 - 4S|)v?2 

+2ATTBiB2f^ + UttB^ip^] (3.1) 

where (p = (p/mpi. This expression contains the Bi at linear and quadratic order - which 
reflects the and H'^ terms rather than a slow roll related truncation - including ^3 and 
higher terms in H does not add terms of the form BiBjB^ to the potential. Secondly, the 
overall scale of the potential is set by Hq, which fixes the energy scale of inflation. Looking 
at the structure of H^, we can see that the first i powers of ip will each have a linear 
contribution from Bi, up to Bjsi. All other contributions involve quadratic combinations 
like BiBj or Bf. In particular, the ip^~^^ through y?^^ terms cannot be fixed independently 
of the lower order terms. 

Near the fiducial point, p <^ 1 and the higher order terms are effectively ignorable, and 
the only relevant contributions to V{(p) arise from the terms that are linear in Bi. Thanks 



to a well-known result, originally due to Lyth |24, |2^, the change of (f) during the course 



*Note that the present day Hubble constant is specified by h, and is not connected with the Ho used 
here. 
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of inflation is a function of r, the tensor:scalar ratio. From our perspective, the two key 
results are 



r = 16e, (3.2) 
A<p = ^V^AN. (3.3) 

The tensor :scalar ratio measures the relative amplitudes of Pji and Ph, while AA^ and A(j) 
measure the number of e-folds that elapse as the field rolls through A(/).^ These relationships 



are derived at lowest order in slow roll and can be extended [25| to higher orders, but for 



the present analysis these expressions are more than adequate. The immediate lesson to 



be drawn from equation 3^ is that if r is well below the current bound, Aip will be less 
than unity even if AA^ is as large as 50. 

Recalling the relationship between e and Bi we can derive 

^ = flAA^. (3.4) 

mpi 47r 

again at lowest order in slow roll. Rewriting H{(p) in terms of the number of efolds that 
elapse since the ko mode crosses the horizon during inflation we find. 



H{AN) ^ Ho 



l + Bi[ -^^ + • • • + Bm+ ' 
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M+l 



(3.5) 



Even with the constraint derived here from WMAPI+All data, we can confidently say that 
e < 0.02 at 0o- If we are only interested in reconstructing of the piece of the potential 
that contributes to cosmologically relevant perturbations, we only need Bi < 1 for this 
expansion to be well behaved. However, if the bound on e is saturated, then increasing 
AA'^ sees H{AN) dominated by higher order terms and the approximation will become 
unreliable. On the other hand, if it turns out that e <C 0.01 then we could imagine pushing 
AA^ up to 60 without having to put tight constraints on the Bi - and in this case the 
end of inflation would presumably be encoded in a divergence of the Bi for large i. This 
analysis is necessarily heuristic, but is allows us to establish a rough estimate for the region 
of validity of this approximation function of e. 

The initial value of e or Bi effectively fixes Acj) during the interval in which observ- 
able modes leave the horizon. For the chains we run here, the constraints on e from the 
WMAPI+All dataset guarantee that A(p < mpi over the 8 e-folds following the moment 
when the ko mode crosses the horizon. However, as observational bounds on r tighten, the 
permissible range of A(p will decrease. Given that B2 is typically very much less than unity, 
the expansion appears to be well behaved in the region of interest, but this expansion is 
almost certainly asymptoically rather than absolutely convergent, and would eventually 
diverge if a sufficiently large number of terms was included. 



^Note that several different definitions for r can be found in the Uterature so one frequently encounters 
this relationship with a different numerical coefficient. 
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4. Markov Chain Monte Carlo Analysis 



We use a Markov Chain Monte Carlo (MCMC) technique ||T|, ||, ||, ||, |4|, |4|] to evaluate 
the likelihood function of model parameters. The MCMC is used to simulate observations 
from the posterior distribution 'P(a|x), of a set of parameters a given event x, obtained 
via Bayes' Theorem, 



Via\x) 



V{x\a)V{a) 
J V{x\a)V{a)da' 



(4.1) 



where V{x\q) is the likelihood of event x given the model parameters a and V{a) is the 
prior probability density. The MCMC generates random draws (i.e. simulations) from the 
posterior distribution that are a "fair" sample of the likelihood surface. From this sample, 
we can estimate all of the quantities of interest about the posterior distribution (mean, 
variance, confidence levels). A properly derived and implemented MCMC draws from the 
joint posterior density V{a.\x) once it has converged to the stationary distribution. 

We use modified versions of the CosmoMC^/CAMB'' packages for our computations, 
utilizing eight chains per model and a conservative convergence criterion [^. For our 
application, a denotes a set of cosmological parameters, the number of which varies with 
the model chosen to describe the primordial power spectrum (see below). We always vary 
4 "late-time" parameters: the physical energy density in baryons Qbh'^, the total physical 
energy density in matter QcDMh'^, the optical depth to reionization r, the Hubble constant 
h in units of 100 km s~^ Mpc~^. The "event" x is the likelihood of obtaining the total 
observational dataset (either WMAPI+All, WMAPII or WMAPII+SDSS in this paper) if 
we assume that the actual universe has an underlying set of parameter values given by a. 

In order to understand the distinctions between the Hubble slow roll parameters and 
the characterization of the spectrum in terms of an index and running, we perform three 
distinct analyses. Firstly, we present results for the traditional power law primordial ob- 
servables: the scalar spectral index, n^, its running with scale, dus/dlnk, the tensor to 
scalar ratio, r, the tensor spectral index, nt, and the scalar spectral index, Ag (all defined 
at the fiducial scale ko = 0.002 Mpc~^) We are including the amplitude and indices for 
both the scalar and tensor spectra. In single field slow roll inflation these values are re- 
lated. However, we wish to avoid the use of any inflationary priors in these chains, but 
many analyses adopt this ansatz in practice since the tensors are currently unobserved. So 
far as we know, this particular parameter set has not be used before in this context - as 
one would expect, parameters corresponding to the (currently unobserved) tensor power 
spectra are very weakly constrained. We denote this formulation by {ns,nt}. 

Secondly, we run chains where the HSR parameters are varied inside of the Markov 
chains, but the spectra are computed by converting these values into "traditional" power 



http : // cosmologist . info/ cosmomc/ 
'^http : // camb . info/ 



law variables, using the formulae below: 

n, = 1 + 2r/ - 4e - 2(1 + C)e^ " ^(3 - 5C)er/ + ^(3 - C)i (4.2) 

r = 16e[l + 2C(e -??)] (4.3) 

rif = -2e-(3 + C)e2 + (l + C)ery (4.4) 
dus 1 \ 2— - 4 



dink l-ey dN dN 

4a.C,.A_>,3-5C,(.^.,A)^,3-C,^} ,4.5, 

where C = 4(ln 2 + 7) — 5, and we note that C differs from the C in the expressions for the 
spectra. We refer to this formulation as ns{e,rj). 

The scale dependence of the slow roll parameters is obtained directly from the flow 
equations, 

^ = M,-^) (4.6) 
II = + { (4.7) 
§^^ari- 2e) + 'A« (4.8) 

which we have truncated at third order by setting ^X^ = 0. In practice we will set 
^ = ^Xh = as well, so that only t] and e vary as the universe evolves. Finally, denotes 
the number of e-folds of inflation, which does not not explicitly appear in this formulation. 

The initial values of the HSR parameters are assumed to be set at the fiducial scale 
ko = 0.002 Mpc~^. The third and final parameterization uses the HSR parameters directly, 
calculating the power spectrum using the formulae 2.12 and p. 13 . The power spectrum is 
normalized at ko by setting 

where eo and rjQ are the values at ko- 



5. Results 

We now turn to the results of our experimentation with different characterizations of the 
primordial spectrum. In all the chains discussed below, we work with the following set of 





Dataset 


Variables 


Comments 


1. 


WMAPI+All 


{ns,nt, ■■■} 


Standard spectral variables 


2. 


WMAPI+All 


{ns{e,r]),nt{e,r]),---} 


Slow roll parameters in chains 


3. 


WMAPI+All 


HSR {e,r]} 


2D HSR 


4. 


WMAPII 


HSR {e,r?} 


2D HSR 


5. 


WMAPII + SDSS 


HSR {e,r?} 


2D HSR 



Table 1: MCMC Parameter fits described in this paper. 
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"late time" parameters, il}yh? , QcDuh^: h and r, where these variables have their usual 
meanings. Without careful tuning, the inflationary universe is extremely close to flat, so it 
is consistent with our imposition of an inflationary prior to set f^Totai = 1- The contribution 
from dark energy is thus Vt\ = 1 — ilj ~ ^CDM, and we take its equation of state to be 
w\ = — 1. We assume instantaneous reionization at an optical depth of r, the HST Key 
Project flat prior on the Hubble parameter |3^, and a flat prior on the age of the universe, 
as described in the Introduction. 

In the ID plots that follow, we display the probability distributions for each parameter 
marginalized over all the other parameters. In the 2D figures, the solid contours are the 
marginalized 1 and 2 a joint probability contours for each pair of parameters. Each plot 
contains data drawn from 8 chains, which were well-converged according to the Gelman- 
Rubin criterion. The CMB spectra were computed with an appropriately modified version 
of CAMB, and the chains themselves run with the CosmoMC package. 

We ran Monte Carlo Markov Chains over the WMAPI+All dataset for a three different 
parametrizations of the fundamental spectral variables, as outlined in Table 1. In Model 
1, the primordial parameters {ng, dns/dlnk,r,nt, As} measured at ko = 0.002 Mpc~^ were 
varied within the chains. This covers much the same ground as many previous analyses (e.g. 
||3l|] ) - the one novelty being that we have fitted to both the amplitude and index of the 
tensor spectrum. As might be expected in the absence of an observed tensor contribution 
to the primordial perturbations, these parameters are only very weakly constrained. The 
one constraint we do make is that nt < 0, so that the tensor spectrum is red.^ However, 
such models also typically have a vanishingly small tensor amplitude at CMB scales, so it 
is a very mild assumption. We then ran fits using the two different slow roll schemes - 
that of Leach et al. |8|, ^, |lO|, and the HSR formalism described in this paper. We see that 
up to second order they give very similar results, as we would have expected. Finally, we 
repeat the process for WMAPII on its own, and WMAPII+SDSS, obtaining constraints 
on e and rj via the HSR formalism. 

5.1 WMAPI+All with an Inflationary Prior 

Looking at Figure |l], we see that the with the inclusion of an arbitrary tensor spectrum 
and in the absence of any inflationary prior, the allowed parameter ranges are very broad 
when compared to the values typically extracted from fits to recent datasets. We can 
attribute this to the contribution of the tensor component at low i, which can be very 
large in these fits. Moreover, we see that a small value of r is correlated with a very red 
tensor spectrum. In this case the tensor modes contribute strongly at the longest scales, 
and the normalization has to be reduced in order to obtain the correct normalization for 
the lowest lying Cis. Likewise, we see evidence (at the 2a level) for a substantial running 
in the scalar index. This mirrors the result for a obtained from the WMAPII dataset p8| . 

*The height of the tensor spectrum scales with the energy density of the universe - and this is strictly 
decreasing unless one has either "phantom" energy with w < — 1 |^^, or if the perturbations are laid down 
during a contracting phase, as in the ekpyrotic and pre-big bang |Q scenarios. These models do have 
blue gravitational spectra, but the tensor amplitude at cosmological scales is typically infinitesimal. 
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Figure 1: The results of three different fits to the WMAPI+All dataset. The blue lines are derived 
from the usual scalar and tensor spectral indices (Model 1). The black lines are the results of a fit 
using the ns(e, ry) approach (Model 2), and the red lines correspond to the HSR algorithm described 
in this paper (Model 3). The top plot in each column shows the probability distribution function 
for each of the different "primordial" variables, which the other plots show their joint 68% and 
95% confidence levels. When we impose an inflationary prior, we convert e and rj back into the 
parameter values at ko. 



Conversely, when we express the spectrum in terms of e and r] we see that the available 
region of parameter space is dramatically reduced - particularly in the case of the tensor 
modes, and dng/dlnk, which reflects the imposition of the slow roll prior. In particular, 
many analyses of the primordial perturbations impose a slow roll prior that relates the 
amplitude and slope of the tensor perturbation spectrum, and doing so dramatically reduces 
the range of the (nj, r) plane that is compatible with the data. Likewise, the natural "center 
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of gravity" for the WMAPI+All dataset prefers a negative value of a, the running of the 
scalar spectral index. Provided we use only the first two terms in the slow roll expansion 
(e and rj), this quantity can be large only if the scalar spectral index is far from unity. 
If the central values of a found from these two different parametrization do not change 
in future datasets, but the error ellipses contract, then this would provide an alternative 
"consistency" check on simple forms of slow roll inflation. Any failure of the data to 
satisfy would demonstrate either than the underlying inflationary model has a non-trivial 
potential, or that the perturbations are inconsistent with the assumption that they are 
generated by an inflationary epoch generated by a single, minimally coupled scalar field. 
We fit the WMAPI+All dataset using both Leach et aZ.'s approach (Model 2) and our 
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Figure 3: The allowed values of e and rj derived from the MCMC chains for Models 3, 4 and 5. 
In each case we show the joint 68% and 95% confidence contours. The blue lines are WMAPI+AU 
(Model 3), the red lines are WMAPII (Model 4) and the black lines are WMAPII+SDSS (Model 5) 



HSR algorithm (Model 3). The parameter values obtained for these fits are shown on their 
own in Figure ^. The two approaches lead to very similar results for the "early universe" 
cosmological parameters. In Figure ^ we show the constraints on e and r/ for the HSR fit. 
The most obvious feature here is the degeneracy direction following from — 1 ~ 2r/ — 4e. 



Parameter 


WMAPH 


WMAPH±SDSS 


ris 


0.982 ±0.020 


0.973 ±0.017 


dus/dln k 


n 002+° °°^ 


U.UUUU_o.0005 


ln[10i%] 


3.07 ±0.09 


3.17 ±0.06 


r 


< 0.465 (95% CL) 


< 0.200 (95% CL) 


nt 


> -0.059 (95% CL) 


> -0.025 (95% CL) 



Table 2: Constraints on primordial parameters, defined at fc = 0.002 Mpc ^. Constraints are at 
the 68% level unless otherwise noted. 
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Figure 4: We show the resuhs for using the HSR formahsm to apply an inflationary prior to 
WMAPll (red) and WMAPII+SDSS (black). We again show joint 68% and 95% confidence con- 
tours. 



5.2 WMAPII and Inflation 

We use the HSR formalism to apply an inflationary prior to the WMAPII (Model 4) 
and WMAPII+SDSS (Model 5) datasets. The constraint on e and r] is shown in Figure 
m, and we immediately see that WMAPII on its own constrains the data slightly more 
strongly than the previous "global" dataset that comprises WMAPI+All. However, the 
more interesting result here is WMAPII+SDSS, which we see leads to a considerably tighter 
set of constraints, and starts to break the degeneracy between e and t] seen in all the other 
datasets. 

We plot the corresponding parameter values in Figure ^, and present the corresponding 
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numerical bounds in Table ^ We immediately see that there are strong constraints on the 

tensor contribution and the running of the scalar index for the WMAPII+SDSS case. 

In both cases these arise via the breaking of the degeneracy that allows e to run out to 

relatively large values. The joint 95% confidence interval for the value of r is lower than is 

obtained for a generic (n^, r) fit on the same dataset |3^. Moreover, the central value of 

has moved slightly towards unity, relative to ||3^, though the distributions are consistent 

at less than the la level. Consequently, while the the pure Harrison-Zel'dovich spectrum 

is disfavoured, it is not strongly excluded. Finally, we note that this data continues to 

be a good for fit for m'^cfP' inflation, and is excluded at more than 2cr-joint by the 

WMAPII+SDSS fit on account of its primordial tensor spectrum alone. 

Note that for the WMAPII data, we use an updated beam error module which differs 

from the version that was used in [ 38 1 . This is available as an option in the public WMAPII 

likelihood software,^ and uses the full off-diagonal Cj^ covariance matrix, a Gaussian+ln 

normal approximation for the Cf^ likelihood, and a fixed fiducial cj^^°^^ to propagate 

the beam errors in the likelihood function. The default options in the likelihood software 

assume a diagonal off-diagonal Cf^ covariance matrix, a Gaussian approximation for the 
jti 



CT^ likelihood, and cj^'^°^^ varied in the steps in the MCMC to propagate the beam error. 



Note that, unlike in [38|, we do not marginalize over the amplitude of the SZ fluctuations. 
This should be a mild assumption since there is no detection of SZ fluctuations, and the 
expected amplitude of SZ contamination in the data is at the level of the weak lensing of 
the CMB which was also neglected. Also, adding such a parameter could cause spurious 
degeneracies within our parameter set, a consequence which we wanted to avoid. 

We show our constraints from WMAPII and the closest corresponding constraints from 
j38U (a fit using r and Ug only, with no running) in Figure ^. The differences in the likelihood 
approximation between our analysis and that of Spergel et al. outlined above only shifts the 
constraints by considerably less than la-marginalized if the (r, n^) parametrization is used 
with the differing beam treatments and the SZ marginalization. Therefore the differences 
in the distributions for and r seen in Fig 5 arise from the imposition of the slow roll 
prior. 



6. Discussion 



The results in the previous section display the effectiveness of the flow reconstruction 
algorithm, and clarify the distinctions between this approach and other methods that have 
been tried in the past. To recapitulate, the primary advantages of this approach is that 
is does not require one to Taylor expand the perturbation spectrum about an arbitrary 
fiducial point, and working with the full inflationary dynamics is guaranteed to capture any 
consistency relationships that relate the spectral parameters. In particular, this process 
represents a step forward from Monte Carlo reconstruction, in that it avoids the use of 
an arbitrary "measure" on the parameter space occupied by the initial values of the flow 
parameters. 

® However, it is turned off by default. 
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Figure 5: Comparison of constraints in the r vs. Us plane: The red hnes correspond to the HSR 
algorithm for WMAPII described in this paper (Model 3) converting e and 77 back into the plotted 
parameter values at fco- The black lines are for a fit to the power law variables r and ris for the 
same data set. In both cases, joint 68% and 95% confidence levels are plotted. 

As stated in the Introduction, a similar approach to the problem of determining the 
values of the slow roll parameters (and from them, the inflationary potential) has been 
proposed by Leach and collaborators ||8|, Their method is based on Taylor expanding 

the slow roll parameters about their values at some fiducial point. For many purposes this 
method is equivalent to the flow equation approach developed here, but it differs in two 
crucial ways. Firstly, the flow equations allow reconstruction to be carried out to arbitrary 
order in the slow roll parameters, and secondly when we use the flow equations we can 
explicitly check that our choice of slow roll parameters is consistent with obtaining the 
required amount of inflation. In particular, unless e is vanishingly small, we typically 
require at least 50 efolds of inflationary growth to occur after the mode /cq = 0.002Mpc~^ 
leaves the horizon. Looking at Figure we see that imposing this requirement actually 
removes a piece of the parameter space that overlaps with the allowed regions for e and r] 
that one obtains from the MCMC fits. We caution against simply imposing an additional 
constraint on the parameter space by insisting that > 55 - however, we can say that 
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0.05 



Figure 6: We plot the difference between the scalar power spectrum computed using ns{r],e) 
and the full HSR parameters at 1 Mpc"'^. The contours are 1% apart, and a positive difference 
corresponds to the HSR result being larger than the ns{e,r]) value. We have overlapped the joint 
68% and 95% confidence intervals for e and rj obtained for WMAPH (red) and WMAPH+SDSS 
(black). The white region corresponds to parameter values which support less than 15 e- folds of 
inflation after the ko mode leaves the horizon. 



models in this apparently excluded region must have a significant contribution from higher 
order terms in the slow roll expansion in order to be viable. 

More generally, we have chosen to expand H{(j)) as a polynomial in (p. While this is 
perhaps a natural choice it is certainly not unique - in principal one could use any set of 
non-degenerate functions of as a basis for the MCMC fit. In particular, if one wanted 
to generalize this method to potentials with "features" it may be useful to introduce a 
different parameterization of H{(f))}^ In this case we may need to replace Eqs. p. 12 and 



"in this context see also the analysis of 
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Figure 7: We show the number of e-folds that elapses between the slow roll equations becoming 
singular (marking the end of inflation) and the mode fco leaving the horizon, along with the WMAPII 
and WMAPII+SDSS constraints, with the same color coding as Figure |^. 



2.15 with a direct integration of the differential equations describing the evolution of the 
scalar and tensor perturbations to calculate their power spectra as a function of k, before 
feeding them into the CMB power spectrum calculation. 

In Figures ^ through |^ we explore our different constraints on e and i].^^ Figure ^ 
shows the difference in the primordial power spectra (not the C^'s themseves) computed 
directly from the HSR expansion, and the values obtained using the ns{€,rj) approach of 
Leach et al. We see that these overlap well, except for relatively extreme values of the 
slow roll parameters. However, as we discuss in the Introduction, the use of heterogeneous 
datasets may eventually lead to constraints on the primordial power spectrum for modes 



^^In all plots, we are showing eo and ?7o, the instantaneous values of the slow roll parameters when k — fco, 
but we drop the subscripts for clarity. 
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Figure 8: dris/dhik is sliown as a function of e and 77, along witli our constraints on e and rj. 

with k considerably smaller than 1 Mpc^^, and in this case the HSR formalism will be 
indispensable. We plan to address this issue in more detail in future work. 

In Figure we show the number of e-folds that elapse after the mode ko leaves the 
horizon, as a function of e and tj.^"^ This result has to be interpreted with care, since while 
higher order terms in the slow roll expansion may have a negligible impact on the spectrum 
at CMB scales, they can exert considerable leverage between the generation of cosmological 
perturbations and the end of inflation. However, this shows it is straightforward to add 
constraints on the total number of e-folds of inflation directly to the parameter estimation 
process. Moreover, we see that a naive imposition of the requirement that > 55 would 

Strictly speaking, we plot the number of e-folds that elapses before the slow roll hierarchy becomes 
singular, with e — > 00. In practice this is effectively synonymous with the end of inflation. The singularity 
has no physical significance, since the HSR expansion breaks down when infiation ends. 
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Figure 9: The inflationary parameter splace is shown as a function of e and 77, along with our 
constraints on e and rj. 



rule out a piece of the parameter space that is otherwise consistent with the WMAPII 
results. Finally, we remind the reader that while it may be significant when a model fails 
to produce the requisite amount of inflation, ii N ^ 55 we cannot rule the model out of 
consideration, as inflation can halt when an instability in the scalar field potential allows 
the field point to roll in an orthogonal direction. 

The running of the spectral index, dus/dlnk, is displayed in Figure ^. This recapitu- 
lates the well known result that the running is typically very small for simple inflationary 
models. If this parameter really is significant (as is tentatively suggested by the WMAPII 
dataset), it will rule out essentially all inflationary models which can be described in terms 
of just two slow roll parameters. At the very least we would need to include higher order 
terms in the slow roll expansion in order to flt such a signal. Further, the large reduction is 
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the allowed values of the running when one adds the SDSS to the fit (as seen in Figure ^ 
is explained here, since the 95% contour for WMAPII corresponds to a much larger value 
(although still tiny) of dng/dlnk than the corresponding contour for WMAPII+SDSS - 
which follows from the quadratic dependence of dug/dlnk on e and rj. Finally, Figure ^ 
shows the division of the inflationary parameter space into large field, small field and hy- 
brid models |53, 54 1, overlaid with our constraints on e and rj. We see that all three classes 
of models are allowed at the la level, but the range of parameter values within each class 
is tightly restricted by the WMAPII+SDSS results. 

It is of particular interest to us that the differences between the constraints found 
from WMAPII and WMAPII+SDSS appear to arise via the dependence of dng/dlnk on 
e and rj. In particular, one can see from Figure |^ that the degeneracy direction in the 
(e, rj) plane defined by constant (at lowest order in slow roll) is broken when one adds 
the running terms. Fortuitously, as one moves out along this line, dus/dlnk > and r 
grows with e. In a conventional fit, we see a degeneracy between the tensor amplitude, 
the scalar amplitude and the scalar spectral index. Increasing r adds power to the Ci in 
the Sachs- Wolfe tail of the spectrum, and Ag is driven down to compensate. However, the 
tensors do not contribute to the temperature anisotropies past the first acoustic peak and 
to large structure, so rises to offset the lower scalar amplitude at small scales. 

By using the HSR equations, we necessarily increase dug/dlnk when we increase e 
and r - and the "lever arm" felt by this term is magnified when we include SDSS data. 
The HSR approach thus implicitly imposes a "higher order" consistency condition, where 
the running of the scalar index is correlated with the tensor amplitude |5|, To our 
knowledge, this is the first time that a cosmological dataset has been detailed enough to 
justify the inclusion of such higher order correlations, underscoring the importance of the 
new WMAPII results. We would add two qualifications to this claim however. The first is 
that we only work with two slow roll parameters. In the usual expansion dng/dln k depends 
on H"'{(j)), so the results may change significantly if we add the ^ term to our parameter 
set.^'^ Secondly, while we know from Spergel et al. [^] that the data is consistent with a 
spectral index that decreases with k, a mildly positive running is ruled out. If the central 
value for the running was closer to zero, the ability of these higher order terms to constrain 
the infiationary parameter space would be reduced. On the other hand, if dus/dlnk is 
close to the central value suggested by [38|, our results would imply that theories based on 
a single, minimally coupled infiaton is inconsistent with observations if the potential can 
be well specified solely in terms of e and rj, ruling out a vast class of simple infiationary 
models. 



7. Conclusions 

In this paper we have introduced a new approach to introducing an inflationary prior into 
a cosmological parameter estimation. We proceed by varying the values of the the Hub- 
ble Slow Roll parameters directly inside the Monte Carlo Markov Chains, and evolving 
the values of these parameters using the "flow equations" of the slow roll hierarchy. This 

^^Or, equivalently, F"'(</>). 
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approach has much in common with Leach et aZ.'s use of the slow roll expansion in cos- 
mological parameter estimations ^, |l^ , and Monte Carlo reconstruction which uses the 
flow hierarchy to generate inflationary models which match a predefined set of constraints 
0, combining the advantages of both techniques. In particular, since we do not need to 
expand the power spectrum or other physical quantities as a Taylor series about some fixed 
point, we expect that this method will be particularly suited to high quality, heterogeneous 
datasets that we expect will become available in the future. 

We can see several directions in which to extend this work. The weak but persistent 
evidence for a large running in the scalar spectrum found in the WMAPII data provides 
a compelling motivation for extending our results to include the third order term in 
the slow roll expansion. From a numerical perspective, this makes the likelihood surface 
significantly more complicated and increases the time taken for the MCMC chains to con- 
verge, but we expect that the problem will be computationally tractable. Moreover, as 
further cosmological datasets emerge, it will obviously be of considerable interest to see 
how the allowed region of HSR parameter space evolves and, hopefully, contracts. From 
a theoretical standpoint, we plan to look more carefully at the use of the HSR flow code 
with high quality heterogeneous datasets that incorporate Lyman-a forest or other short 
wavelength probes of the primordial spectrum. In particular, this approach would be very 
well suited to a dataset that contained a direct detection of primordial gravitational waves 
from BBO or a similar experiment, and we suspect this would allow us to put exquisite 
constraints on the form of the inflationary potential. 

On a practical level, we have used the Hubble Slow Roll equations to apply an inflation- 
ary prior to both the WMAPI+All dataset, a compendium of several CMB experiments, 
and the new WMAPII data, with and without the SDSS power spectrum. In doing so, we 
find that the central value of Ug moves slightly towards the scale free Harrison-Zel'dovich 
result. The discrepancy between this and the results in |Q is explained by the use of 
slightly different assumptions in our likelihood function, and sheds light on the possible 
systematic errors that arise when fitting cosmological parameters to this accuracy. There 
are other small systematic effects which can move the centroid of the distribution to- 
wards redder values, and in particular, a stronger upper limit on r (e.g. by using a larger 
complement KaQVW of WMAP bands in the polarization analysis) will strengthen the 
case for deviation from scale invariance by cutting down the residual ns,T degeneracy, as 
will dropping the assumption of instantaneous reionization (David Spergel, private com- 
munication). This is clearly a question which more data will settle. On the other hand, we 
find a tighter constraint on the tensor amplitude, r, than provided by [p8| . 

All of the major classes of single field inflationary models are allowed by our con- 
straints, although many speciflc models are ruled out. However, we find evidence that the 
combination WMAPII+SDSS is powerful enough to be sensitive to terms that are second 
order in the slow roll parameters e and rj, via the contribution of these terms to the run- 
ning of the scalar spectral index, dus/dlnk. This evidence is tentative and could easily 
evaporate as the overall cosmological dataset continues to improve, but it is indicative of 
the usefulness of the Hubble slow roll formalism as a tool for analyzing cosmological data. 
Finally, and independently of our detailed results, we can conclude that the era of preci- 
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sion cosmology is now in full flight, and that the promise that astrophysical data will put 
meaningful constraints on the inflationary potential is beginning to be fulfilled. 
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